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ABSTRACT 

Q . We present a detailed analysis of 101 white dwarf-main sequence binaries (WDMS) 

■ from the Sloan Digital Sky Survey (SDSS) for which multiple SDSS spectra are avail- 
C/3 ' able. We detect significant radial velocity variations in 18 WDMS, identifying them 

as post common envelope binaries (PCEBs) or strong PCEB candidates. Strict upper 
limits to the orbital periods are calculated, ranging from 0.43 to 7880 d. Given the 
sparse temporal sampling and relatively low spectral resolution of the SDSS spectra, 
our results imply a PCEB fraction of >15 % among the WDMS in the SDSS data 
f — ' base. Using a spectral decomposition/fitting technique we determined the white dwarf 

. effective temperatures and surface gravities, masses, and secondary star spectral types 

for all WDMS in our sample. Two independent distance estimates are obtained from 
the flux scaling factors between the WDMS spectra, and the white dwarf models and 
. main sequence star templates, respectively. Approximately one third of the systems 

in our sample show a significant discrepancy between the two distance estimates. In 
t^- , the majority of discrepant cases, the distance estimate based on the secondary star is 

■ too large. A possible explanation for this behaviour is that the secondary star spec- 
tral types that we determined from the SDSS spectra are systematically too early by 
1-2 spectral classes. This behaviour could be explained by stellar activity, if cover- 

■ ing a significant fraction of the star by cool dark spots will raise the temperature of 

the inter-spot regions. Finally, we discuss the selection effects of the WDMS sample 
provided by the SDSS project. 

Key words: accretion, accretion discs - binaries: close - novae, cataclysmic variables 
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1 INTRODUCTION 

A large fraction of all stars in the sky are part of binary 
or multiple systems (Ibcn 1991). If the initial separation of 
the main-sequence binary is small enough, the more mas- 
sive star will engulf its companion while evolving into a 
red giant, and the system enters a commo n envelope (CE, 
e.g. lLivio fc Sokerll 19881 ; llben fc Livioll 19951 ). Friction within 
the CE leads to a rapid decrease of the binary separation 
and orbital period, and the energy and angular momentum 
extracted from the binary orbit eventually ejects the CE. 
Products of CE evolution include a wide range of important 
astronomical objects, such as e.g. high- and low-mass X-ray 
binaries, double degenerate white dwarf and neutron star bi- 
naries, cataclysmic variables and super-soft X-ray sources - 
with some of those objects evolving at later stages into type 



la supernova and short gamma-ray bursts. While the con- 
cept of CE evolution is simple, its details are poorly under- 
stood, and are typically described by p a rametrised models 
llPaczvnskil Il97rj ; iNelemans et all l200d ; iNelemans fc ToutJ 
l2005f ). Consequently, population models of all types of CE 
products are subject to substantial uncertainties. 

Real progress in our understanding of close binary evo- 
lution is most likely to arise from the analysis of post com- 
mon envelope binaries (PCEBs) that are both numerous and 
well-understood in terms of their stellar components - such 
as PCEBs containing a white dwarf and a main sequence 
staiQ. Wh ile detailed population models are already avail- 
able, (e.g. lWillems fc Kolblliooi ). there is a clear lack of ob- 

1 Throughout this paper, we will use the term WDMS to refer 
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Figure 1. The spectrum of SDSS J005208.42-005134.6, a typi- 
cal WDMS in the SDSS data base. The white dwarf is clearly 
visible in the blue while the low mass companion dominates the 
red part of the spectrum. Evident are the Ho emission line, and 
the Nal AA 8183.27,8194.81 absorption doublet, originating on the 
companion star. These features are shown in the small insets on 
a velocity scale, and are used to measure the radial velocities of 
101 WDMS for which multiple SDSS spectra exist in DR5. See 
also Fig. [2] 

servational constraints. ISchreiber fc Gansickd (2003) showed 
that the sample of well-studied PCEBs is not only small, 
but being drawn mainly from "blue" quasar surveys, it is 
also heavily biased towards young systems with low-mass 
secondary stars - clearly not representative of the intrinsic 
PCEB population. 

The Sloan Digital Sky Survey (SDSS) is cur- 
rently providing the possibility of dramatically improv- 
ing the observational side of PCEB studies, as it has 
already identified close to 1000 WDMS (see Fig.QJ 
with hundreds more to follow in futu r e data releases 
jRavmond et al.ll2003l; ISilvestri et al.ll2006l ; lEisenstein et all 
l200d : ISouthworth et alJl2007t ). Within SEGUE, a dedicated 
program to identify W DMS containing cold white dwarfs is 
successfully underway (|Schreiber et al.ll2007l ). 

Identifying all PCEBs among the SDSS WDMS, and 
determining their binary parameters is a significant obser- 
vational challenge. Here, we make use of SDSS spectroscopic 
repeat observations to identify 18 PCEBs and PCEB can- 
didates from radial velocity variations, which are excellent 
systems for in-depth follow-up studies. The structure of the 
paper is as follows: we describe our WDMS sample and the 
methods used to determine radial velocities in Sect .[2] In 
Sect. [3] we determine the stellar parameters of the WDMS 
in our sample. In Sect .[4] we discuss the fraction of PCEBs 
found, the distrib ution of stellar param eters , compare our 
result s to those of iRavmond et all (|2003l ) and lSilvestri et all 
(2006), discuss the incidence of stellar activity on the sec- 
ondary stars in WDMS, and outline the selection effects of 
SDSS regarding WDMS with different types of stellar com- 
ponents. 



to the total class of white dwarf plus main sequence binaries, and 
PCEBs to those WDMS that underwent a CE phase. 
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Figure 2. Fits to the Nal AA 8183.27,8194.81 absorption doublet 
(right panels) and the Ha emission line (left panels) in the four 
SDSS spectra of the WDMS SDSS J024642. 55+004137.2. The 
SDSS spectroscopic identifiers (MJD, Plate-ID and Fibre- ID) are 
given in the top left corner of the Ho panels. Na I has been fit- 
ted with a double-Gaussian of fixed separation plus a parabola, 
Ho with a Gaussian plus a parabola. In this system, radial ve- 
locity variations are already obvious to the eye. The top three 
spectra are taken in a single night, the bottom one is combined 
from data taken on three nights, MJD = 52970, 52972, and 52973. 
The widths of the Gaussians fitting the Nal doublet are (top to 
bottom) 4.6 A, 5.8 A, 5.3 A, and 6.0 A. 



2 IDENTIFYING PCEBS IN SDSS 

SDSS operates a custom-built 2.5 m telescope at Apache 
Point Observatory, New Mexico, to obtain ugriz imag- 
ing with a 120-megapixel camera covering 1.5 deg 2 at 
once. Based on colours and morphology, objects are then 
flagged for spectroscopic follow-up using a fibre-fed spectro- 
graph. Each "spectral plate" refers physically to a metal 
plate with holes drilled at the positions of 640 spectro- 
scopic plus calibration targets, cov ering ~ 7d e g 2 . Te chni- 
cal details on S DSS are given by lYork et alj <|2000h and 
IStoughton et al l (|2002l ). The main aim of SDSS is the iden- 
tiflc ation of galaxies (e.g. IStrauss et all 12002] ) and quasars 
(e.g. lAdelman-McCarthv et aLlbOOd lTwith a small number 
of fibres set aside for o ther projects, e.g. find ing cataclysmic 
variables and WDMS i|Ravmond et alfeoOot ). 

A feature of SDSS hitherto unexplored in the study of 
WDMS is the fact that ~ 10 per cent of the spectroscopic 
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SDSS objects are observed more than onc^: the detection 
of radial velocity variations between different SDSS spectra 
of a given WDMS will unambiguously identify such a sys- 
tem as a PCEB, or a strong PCEB candidate. Throughout 
this paper, we define a PCEB as a WDMS with an up- 
per limit to its orbital period < 300 d, a PCEB candidate 
as a WDM S with periods 300 d < P olh < 1500 d, following- 
Fig. 10 from lWillems fc Kolbl l|2004 ). which shows the period 
and mass distribution of the present-day WDMS population 
at the start of the WDMS binary phase. WDMS with pe- 
riod > 1500 d have too large binary separations to undergo a 
CE phase, and remain wide systems. While these definitions 
depend to some extend on the detailed configuration of the 
progenitor main seque nce binary, the population model of 
IWillems fc Kolbl (|2004h predicts a rather clean dichotomy. 

We have searched the DR5 spectroscopi c data base for 
multip le e xposures of all the WDM S listed bv lSilvestri et all 
(|2006f ) and lEisenstein et al.l i|2006l ), as well as a set of WDMS 
independently found in the SDSS data by our team. This 
search resulted in a sample of 130 WDMS with two to seven 
SDSS spectra. Among those WDMS, 101 systems have a 
clearly pronounced Nai AA 8183.27,8194.81 absorption dou- 
blet and/or Ha emission in their SDSS spectrqj, and were 
subjected to radial velocity measurements using one or both 
spectral features. The Nai doublet was fitted with a sec- 
ond order polynomial and double-Gaussian line profile of 
fixed separation. Free parameters were the amplitude and 
the width of each Gaussian and the velocity of the dou- 
blet. Ha was fitted using a second order polynomial plus 
a single Gaussian of free velocity, amplitude and width. We 
computed the total error on the radial velocities by quadrat- 
ically adding the uncertainty in the zeropoint of the SDSS 
wavelength calibration (10 km s -1 , IStoughton et al.l |200"^ ) 
and the error in the position of the Nan/Ha lines deter- 
mined from the Gaussian fits. Figure[2] shows the fits to the 
four SDSS spectra of SDSS J024642. 55+004137.2, a WDMS 
displaying an extremely large radial velocity variation iden- 
tifying it as a definite PCEB. This figure also illustrates an 
issue encountered for a handful of systems, i.e. that the Ha 
and Nai radial velocities do not agree in the latest spec- 
trum (Table[T}. This is probably related to the inhomoge- 
neous distribution of the Ha emission over the surface of 
the companion star, and will be discussed in more detail 
in Sect. 14. ll In total, 18 WDMS show radial velocity varia- 
tions among their SDSS spectra at a 3cr level and qualify as 
PCEBs or strong PCEB candidates. Their radial velocities 
are listed in Table[T]and illustrated in Fig.0and Fig.fi] Three 
systems (SDSS0251-0000, SDSS1737+5403, and SDSS2345- 
0014) are subject to systematic uncertainties in their radial 
velocities due to the rather poor spectroscopic data. The 
radial velocities for the remaining 83 WDMS that did not 
show any significant variation are available in the electronic 
edition of the paper (see Tablef2} • 

We note that special care needs to be taken in estab- 
lishing the date and time when the SDSS spectra where 

2 SDSS occasionally re-observes entire spectral plates, where all 
targets on that plate get an additional spectrum, or has plates 
which overlap to some extent, so that a small subset of targets 
on each plate is observed again. 

3 The SDSS spectra are corrected to heliocentric velocities and 
provided on a vacuum wavelengths scale. 



obtained: a significant fraction of SDSS spectra are com- 
bined from observations taken on different nights (which we 
will call "sub-spectra" in what follows) in which case the 
header keyword MJDLIST will be populated with more than 
one date. The headers of the SDSS data provide the ex- 
posure start and end times in International Atomic Time 
(TAI), and refer to the start of the first spectrum, and the 
end of the last spectrum. Hence, a meaningful time at mid- 
exposure can only be given for those SDSS spectra that were 
obtained in a single contiguous observation. 

A crucial question is obviously how the fact that some of 
the spectra in our sample are actually combinations of data 
from several nights impacts our aim to identify PCEBs via 
radial velocity variations. To answer this question, we first 
consider wide WDMS that did not undergo a CE phase, i.e. 
binaries with orbital periods of > years. For these systems, 
sub-spectra obtained over the course of several days will 
show no significant radial velocity variation, and combining 
them into a single spectrum will make no difference except of 
increasing the total signal-to-noise ratio (S/N). In contrast 
to this, for close binaries with periods of a few hours to a 
few days, sub-spectra taken on different nights will sample 
different orbital phases, and the combined SDSS spectrum 
will be a mean of those phases, weighted by the S/N of the 
individual sub-spectra. In extreme cases, e.g. sampling of the 
opposite quadrature phases, this may lead to smearing of the 
Nai doublet beyond recognition, or end up with a very broad 
Ha line. This may in fact explain the absence / weakness of 
the Nai doublet in a number of WDMS where a strong Nai 
doublet would be expected on the basis of the spectral type 
of the companion. In most cases, however, the combined 
SDSS spectrum will represent an "effective" orbital phase, 
and comparing it to another SDSS spectrum, itself being 
combined or not, still provides a measure of orbital motion. 
We conclude that the main effect of the combined SDSS 
spectra is a decreased sensitivity to radial velocity varia- 
tions due to averaging out some orbital phase information. 
Figure[2] shows an example of a combined spectrum (bottom 
panel), which contains indeed the broadest Nai lines among 
the four spectra of this WDMS. 

In order to check the stability of the SDSS wavelength 
calibration between repeat observations, we selected a to- 
tal of 85 F-type stars from the same spectral plates as our 
WDMS sample, and measured their radial velocities from 
the Can AA 3933.67,3968.47 H and K doublet in an anal- 
ogous fashion to the Nai measurement carried out for the 
WDMS. None of those stars exhibited a significant radial ve- 
locity variation, the maximum variation among all checked 
F-stars had a statistical significance of 1.5er. The mean of the 
radial velocity variations of these check stars was found to be 
14.5 kms" 1 , consistent with the claimed lOkrns -1 accuracy 
of the zero-point of the wavelength calibra tion for the spec- 
tra fr om an individual spectroscopic plate (jStoughton et al] 
2002). In short, this test confirms that the SDSS wavelength 
calibration is stable in time, and, as anticipated above that 
averaging sub-spectra does not introduce any spurious ra- 
dial velocity shifts for sources that have no intrinsic radial 
velocity variation (as the check stars are equally subject to 
the issue of combining exposures from different nights into a 
single SDSS spectrum) . We are hence confident that any sig- 
nificant radial velocity variation observed among the WDMS 
is intrinsic to the system. 
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Figure 3. Radial velocities obtained from the Nal absorption doublet. WDMS J? 3<r RVs variation, i.e. PCEBs, are shown in black. 
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Figure 4. Same as Fig. [3] but for the Ho radial velocities. 



3 STELLAR PARAMETERS 

The spectroscopic data provided by the SDSS project are 
of sufficient quality to estimate the stellar parameters of 
the WDMS presented in this paper. For this purpose, we 
have developed a procedure which decomposes the WDMS 
spectrum into its white dwarf and main sequence star com- 
ponents, determines the spectral type of the companion by 
means of template fitting, and derives the white dwarf ef- 
fective temperature (T e s) and surface gravity (logg) from 
spectral model fitting. Assuming an empirical spectral type- 
radius relation for the secondary star and a mass-radius 
relation for the white dwarf, two independent distance es- 
timates are calculated from the flux scaling factors of the 
template/model spectra. 

In the following sections, we describe in more detail the 
spectral templates and models used in the decomposition 
and fitting, the method adopted to fit the white dwarf spec- 
trum, our empirical spectral type-radius relation for the sec- 



ondary stars, and the distance estimates derived from the 
fits. 



3.1 Spectral templates and models 

In the course of decomposing/fitting the WDMS observa- 
tions, we make use of a grid of observed M-dwarf templates, 
a grid of observed white dwarf templates, and a grid of 
white dwarf model spectra. High S/N ratio M-dwarf tem- 
plates matching the spectral coverage and resolution of the 
WDMS data were produced from a few hundred late-type 
SDSS spectra from DR4. Th ese spectra were c l assifie d us- 
ing the M-dwarf templates of iBeuermann et al.l (|l998l ). We 
averaged the 10 — 20 best exposed spectra per spectral sub- 
type. Finally, the spectra were scaled in flux to match the 
surface brightness at 7500 A and in the TiO abso rption 
band near 7165 A, as defined by IBeuermann! {2006). Re- 
cently, iBochanski et al.l (|2007T l published a library of late 
type stellar templates. A comparison between the two sets 
of M-dwarf templates did not reveal any significant dif- 
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Table 1. Radial velocities of our 18 PCEBs and PCEB candidates, measured from the Ha emission line and/or the Na I AA 8183.27,8194.81 
absorption doublet. The HJDs for SDSS spectra that have been combined from exposures taken in several different nights (see Sect. [2]l 
are set in italics. PCEB candidates with uncertain radial velocity measurements are indicated by colons preceding and trailing the object 
name. Upper limits of the orbital periods are also provided (Sect. l4~2"t . The two spectral components identified in the spectra are coded as 
follows. DA = white dwarf with clearly visible Balmer lines; DC = clearly visible blue continuum without noticeable structure; blx = weak 
blue excess; dM = M-dwarf. 
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-167.2 


± 


11. 


.8 




0246+0041 


2451871.2731 


-95.5 


± 


10.2 


-99.3 


± 


11.1 


2.5 


1724+5620 


2451812.6712 


125.6 


± 


10.2 


160.6 


± 


18. 


.4 


0.43 


DA/dM 


2452177.4531 


163.1 


± 


10.3 


167.2 


± 


11.3 




DA/dM 


2451818.1149 


108.3 


± 


11.1 




± 










2452965.2607 


140.7 


± 


10.8 


135.3 


± 


11.0 






2451997.9806 


-130.6 


± 


10.3 


-185.5 


± 


20. 


1 






2452971.7468 


64.0 


± 


10.5 


125.7 


± 


11.3 




1726+5605 


2451812.6712 


-44.3 


± 


16.7 


-38.9 


± 


12. 


9 


29 


:0251-0000: 


24521744732 


4.1 


± 


33.5 


0.0 


± 


15.4 


0.58 


DA/dM 


2451993.9805 


46.6 


± 


14.6 


47.3 


± 


12. 


5 




DA/dM 


2452177.4530 


-139.3 


± 


24.6 


15.8 


± 


18.3 




: 1737+5403: 


2451816.1187 








-123.5 


± 


28. 


6 


6.6 


0309-0101 


2451931.1241 


44.8 


± 


13.2 


31.5 


± 


13.7 


153 


DA/dM 


2451999.4602 








44.0 


± 


21. 







DA/dM 


2452203.4500 


51.2 


± 


14.1 


48.7 


± 


24.4 




2241+0027 


2453261.2749 


9.1 


± 


17.9 


22.0 


± 


12. 


.4 


7880 




2452235.2865 


27.4 


± 


14.1 


76.3 


± 


16.2 




DA/dM 


2452201.1311 


-60.3 


± 


12.7 


8.1 


± 


12. 


2 






2452250.2457 


28.8 


± 


15.0 


8.1 


± 


33.0 




2339-0020 


2453355.5822 


-29.2 


± 


10.4 


-27.1 


± 


12. 


3 


120 




2452254.2052 


53.9 


± 


11.8 


55.7 


± 


14.3 




DA/dM 


2452525.3539 


-93.6 


± 


12.3 


-90.1 


± 


12. 


.7 






2452258.2194 


15.5 


± 


13.0 


27.9 


± 


19.9 




:2345-0014: 


2452524.3379 








-141.5 


± 


22. 


9 


9.5 




2453383.6493 


50.7 


± 


11.0 


56.5 


± 


12.8 




DA/dM 


2453357.5821 








-19.8 


± 


19. 


3 




0314-0111 


2451931.1242 


-41.6 


± 


12.4 


-51.7 


± 


12.4 


1.1 


2350-0023 


2451788.3516 


-160.3 


± 


16.6 










0.74 


DC/dM 


2452202.3882 


35.6 


± 


10.9 


35.2 


± 


14.4 




blx/dM 


2452523.3410 


154.4 


± 


31.3 













Notes on individual systems. 0246+0041, 0314-0111, 2241+0027, 2339-0020: variable Ha equivalent width (EW); 0251-0000: faint, weak 
Ha emission with uncertain radial velocity measurements; 1737+5403, 2345-0014: very noisy spectrum; See additional notes in Table|5] 

Table 2. 83 WDMS in our sample that did not show a significant a significant radial velocity variation between their different SDSS 
spectra. The complete table is available in the electronic edition of the paper. The first column gives the SDSS object name, the second 
the HJD of the spectrum, in the third column we quote with y and n those spectra which are composed of subspectra taken in different 
nights, the fourth and fith columns provide the Na I AA 8183.27,8194.81 absorption doublet and Ha emission radial velocities, respectively. 
Blanck spaces indicate that no radial velocity measurement could be obtained. 



Object HJD Sub. RV(Ha)km S - 1 RV(Na)kms- 



SDSSJ001247.18+001048.7 


2452518.4219 


y 


0.4 


± 


26.0 


25.6 


± 


13.8 


SDSSJ001247. 18+001048.7 


2452519.3963 


y 


34.1 


± 


46.7 


9.1 


± 


20.0 


SDSSJ001726.63-002451.1 


2452559.2853 


y 


-4.0 


± 


16.2 


-34.3 


± 


10.2 


SDSSJ001726.63-002451.2 


2452518.4219 


y 


-34.1 


± 


11.5 


-32.0 


± 


12.9 


SDSSJ001749.24-000955.3 


2451794.7902 


n 


-47.6 


± 


13.6 


-44.3 


± 


15.7 


SDSSJ001749.24-000955.3 


2452518.4219 


y 


-7.3 


± 


14.9 


-8.7 


± 


11.6 


SDSSJ001855. 19+002134.5 


2451816.3001 


y 


54.2 


± 


17.6 








SDSSJ001855. 19+002134.5 


2451892.5884 


n 


11.0 


± 


16.3 









ference. We also compiled a library of 490 high S/N DA 
white dwarf spectra from DR4 covering the entire observed 
range of T c ff and logp. As white dwarfs are blue objects, 
their spectra suffer more from residual sky lines in the I- 
band. We have smoothed the white dwarf templates at wave- 
lengths > 7000 A with a five-point box car to minimise the 
amount of noise added by the residual sky lines. Finally, we 



computed a grid of synthetic DA white dw arf spectra us- 
ing th e model atmosphere code described by iKoester et al.l 
(2005), covering logg = 5.0 — 9.5 in steps of 0.25 and 
T c ff = 6000 — 100000 K in 37 steps nearly equidistant in 
log(Toff). 
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gf 0.4 



SDSSJ204431. 43-06 1440. 2 
53269- 



-061-140.2 _ 
-1916-166 . 





6000 7000 8000 9000 

Wavelength [A] 



6 0.15 




SDSS.I005457. 61 -002517.0 
51913-394-110 




_j i i , ^ , aL i , , , j 

6000 7000 8000 9000 

Wavelength [A] 



Figure 5. Two-component fits to the SDSS WDMS spectra. Shown are examples for objects with either the M-dwarf or the white 
dwarf dominating the SDSS spectrum. The top panels show the WDMS spectrum as black line, and the two templates, white dwarf and 
M-dwarf, as dotted lines. The bottom panels show the residuals from the fit. The SDSS spectrum identifies MJD, PLT and FIB are given 
in the plots below the object names. 
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Figure 6. Spectral model fits to the white dwarf components of the two WDMS shown in Fig. 5, obtained after subtracting the best-fit 
M-dwarf template. Top left panels: best-fit (black lines) to the normalised H/3 to He (gray lines, top to bottom) line profiles. Top right 
panels: 3, 5, and 10a \ 2 contour plots in the T c ff — logg plane. The black contours refer to the best line profile fit, the red contours to the 
fit of the whole spectrum. The dashed line indicates the occurrence of maximum H/3 equivalent width. The best "hot" and "cold" line 
profile solutions are indicated by black dots, the best fit to the whole spectrum is indicated by a red dot. Bottom panels: the residual 
white dwarf spectra resulting from the spectral decomposition and their flux errors (gray lines) along with the best-fit white dwarf model 
(black line) to the 3850-7150 A wavelength range (top) and the residuals of the fit (gray line, bottom). The T e g and log g values listed in 
Table[5]are determined from the best line profile fit. The fit to the whole spectrum is only used to select between the "hot" and "cold" 
line fit. 



3.2 Spectral decomposition and typing of the 
secondary star 

Our approach is a two-step procedure. In a first step, we 
fitted the WDMS spectra with a two-component model and 
determined the spectral type of the M-dwarf. Subsequently, 
we subtracted the best-fit M-dwarf, and fitted the resid- 
ual white dwarf spectrum (Sect. 133)1 , We used an evolution 



strategy (jRechenberdl 1994h to decompose the WDMS spec- 
tra into their two individual stellar components. In brief, this 
method optimises a fitness function, in this case a weighted 
X 2 , and allows an easy implementation of additional con- 
straints. Initially, we used the white dwarf model spectra and 
the M-dwarf templates as spectral grids. However, it turned 
out that the flux calibration of the SDSS spectra is least 
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reliable near the blue end of the spectra, and correspond- 
ingly, in a number of cases the \ 2 °f the two-component 
fit was dominated by the poor match of the white dwarf 
model to the observed data at short wavelengths. As we 
are in this first step not yet interested in the detailed pa- 
rameters of the white dwarf, but want to achieve the best 
possible fit of the M-dwarf, we decided to replace the white 
dwarf models by observed white dwarf templates. The large 
set of observed white dwarf templates, which are subject to 
the same observational issues as the WDMS spectra, pro- 
vided in practically all cases a better match in the blue part 
of the WDMS spectrum. From the converged white dwarf 
plus dM template fit to each WDMS spectrum (see Fig. 5), 
we recorded the spectral type of the secondary star, as well 
as the flux scaling factor between the M-star template and 
the observed spectrum. The typical uncertainty in the spec- 
tral type of the secondary star is ±0.5 spectral class. The 
spectral types determined from the composite fits to each 
individual spectrum are listed in Table[5]for the PCEBs in 
the analysed sample, and in the electronic edition of this 
paper for the remaining WDMS (see Table|3|. Inspection of 
those tables shows that for the vast majority of systems, the 
fits to the individual spectra give consistent parameters. We 
restricted the white dwarf fits to WDMS containing a DA 
primary, consequently no white dwarf parameters are pro- 
vided for those WDMS containing DB or DC white dwarfs. 

3.3 White dwarf parameters 

Once the best-fit M-dwarf template has been determined 
and scaled appropriately in flux, it is subtracted from the 
WDMS spectrum. The residual white dwarf spectrum is 
then fitted with the grid of DA models described in Sect. l3.ll 
Because of the uncertainties in the flux calibration of the 
SDSS spectra and the flux residuals from the M-star sub- 
traction, we decided to fit the normalised H/3 to He lines 
and omitted Ha where the residual contamination from 
the secondary star was largest. While the sensitivity to the 
surface gravity in creases for the higher Balmer lines (e.g. 
iKepler et aT]|2006h . we decided not to include them in the 
fit because of the deteriorating S/N and the unreliable flux 
calibration at the blue end. We determined the best-fit T e ft 
and log g from a bicubic spline interpolation to the \ 2 values 
on the T c ff — log g grid defined by our set of model spectra. 
The associated la errors were determined from projecting 
the contour at A\ 2 = 1 with respect to the \ 2 of the best 
fit onto the T e g and logg axes and averaging the resulting 
parameter range into a symmetric error bar. 

The equivalent widths of the Balmer lines go through a 
maximum near T a g — 13 000 K, with the exact value being a 
function of \ogg. Therefore, T e fr and logg determined from 
Balmer line profile fits are subject to an ambiguity, often 
referred to as "hot" and "cold" solutions, i.e. fits of similar 
quality can be achieved on either side of the temperature 
at which the maximum equivalent width is occurring. We 
measured the H/3 equivalent width in all the model spectra 
within our grid, and fitted the dependence of the tempera- 
ture at which the maximum equivalent width of H/3 occurs 
by a second-order polynomial, 



Tcff(EW[H/3] ; 



20361 - 3997 log g + 390(log gf 



(1) 



Parallel to the fits to the normalised line profiles, we iBerger et alj 



fit the grid of model spectra to the white dwarf spectrum 
over the wavelength range 3850 - 7150 A (see Fig. |5). The 
red end of the SDSS spectra, where the distortion from the 
M-dwarf subtraction is strongest is excluded from the fit. 
We then use the T e ft and loggr from the fits to the whole 
spectrum, continuum plus lines, to select the "hot" or "cold" 
solution from the line profile fits. In the majority of cases, 
the solution preferred by the fit to the whole spectrum has a 
substantially lower x 2 than the other solution, corroborating 
that it is likely to be the physically correct choice. In a few 
cases, the best-fit T e fi and loggr from the whole spectrum 
are close to the maximum equivalent width given by Eq[T] 
so that the choice between the two line profile solutions is 
less well constrained. However, in most of those cases, the 
two solutions from the line profile fits overlap within their 
error bars, so that the final choice of T c g and log g is not too 
badly affected. 

Once that T e g and loggr are determin ed from the bes t 
line profile fit, we use an updated version of lBergeron et all 's 
(1995) tables to calculate the mass and the radius of the 
white dwarf. Tabled reports T e g, logg, and the white dwarf 
masses for the PCEBs in our sample, while the results for 
the remaining WDMS can be found in the electronic edi- 
tion (TableUJl. We have carefully inspected each individual 
composite fit, and each subsequent fit to the residual white 
dwarf spectrum, and are confident that we have selected the 
correct solution in the majority of cases. Some doubt re- 
mains primarily for a few spectra of very low signal-to-noise 
ratio. The fact that we have analysed at least two SDSS 
spectra for each system allows us to assess the robustness 
of our spectral decomposition/fitting method. Inspection of 
Table[5] shows that the system parameter of a given system, 
as determined from several different SDSS spectra, gener- 
ally agree well within the quoted errors, confirming that our 
error estimate is realistic. 



3.4 An empirical spectral type-radius relation for 
M stars 

In order to use the flux scaling factor between the ob- 
served WDMS spectra and the best-fit M-dwarf tem- 
plates for an estimate of the distance to the system 
(Sect. [33}, it is necessary to assume a radius for the 
secondary star. Since we have determined the spectral 
types of the companion stars from the SDSS spectra 
(Sect. [321)) we require a spectral type-radius relation (Sp — 
R) for M-dwarfs. The community working on cataclysmic 
variables has previously h a d interest in such a rela - 
tion (e.g. iMateo et al.lll985l; ICaillault fc Patterson] Il990h . 
but while Baraffe fc Chabrier ( 19961 ) derived theoretical 
mass/radius/effective temperature-spectral type relation- 
ships for single M-dwarfs, relatively little observational work 
along these lines has been carried out for field low mass 
stars. In contrast to this, the number of low mass stars 
with accurate mass and radius measurements has signifi- 
ca ntly increase d over the past few years (see e.g. the review 
by iRibasI 120061 ), and it appears that for masses below the 
ful ly convective bounda ry stars follow the theoretical models 
bv lBaraffe et al.1 l|l998t ) relatively well. However, for masses 
> 0.3 Mq, o bserved radii exceed the predicted ones. Stellar 
Lopez-M orales 20071) or metallicity effects (e.g. 



activity (e.g. 



20061 ) were identified as possible causes. 
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Table 4. Stellar parameters of the remaining WDMS identified in our sample, as determined from spectral modelling. The complete table 
can be found in the electronic edition. Given are, from left to right: SDSS object name, MJD, fiber and plate number of the spectrum, 
white dwarf effective temperature and error, white dwarf surface gravity and error, mass of the white dwarf and error, distance to the 
white dwarf and error, spectral type of the sec ondary star, distance to th e secondary and error, flag (we refer by and s and e those 
systems which have been studied previously by ISilvestri et al. (2006) and lEisenstein et al" by re those systems whose binary 

components are resolved), and notes. 



Object 


MJD 


plate 


fiber 


T(k) 


err 


log 9 


err 


M(Mq) 


err 


dwd(pc) 


err 


Sp 


dsec(pc) 


err 


flag notes 


SDSSJ000442.00-002011.6 


51791 
52943 


387 
1539 


24 
21 






















2187 
2330 


236 
251 


re 1 


SDSSJ001029.87+003126.2 


51793 


388 


545 


















2 


1639 


339 


s,re 




52518 


687 


347 


13904 


3751 


8.43 


1.16 


0.88 


0.62 


781 


605 


2 


1521 


314 




SDSSJ001247.18+001048.7 


52518 


687 


395 


18542 


5645 


8.75 


0.80 


1.07 


0.41 


661 


446 


3 


830 


132 


e 




52519 


686 


624 


32972 


7780 


8.61 


1.11 


1.01 


0.54 


1098 


930 


3 


936 


119 




SDSSJ001749.24-000955.3 


51795 


389 


112 


72136 


3577 


8.07 


0.14 


0.77 


0.07 


532 


60 


2 


684 


112 


s,e 




52518 


687 


109 


69687 


4340 


7.61 


0.20 


0.59 


0.07 


784 


127 


2 


659 


136 




SDSSJ001726.63-002451.1 


52559 


1118 


280 


12828 


2564 


8.00 


0.46 


0.61 


0.29 


422 


120 


4 


579 


172 


s,e 




52518 


687 


153 


13588 


1767 


8.11 


0.38 


0.68 


0.24 


424 


106 


4 


522 


155 




SDSSJ001855.19+002134.5 


51816 


390 


385 


















3 


1186 


189 






51900 


390 


381 


14899 


9266 


9.12 


1.03 


1.26 


0.54 


445 


330 


3 


1249 


199 






52518 


687 


556 


10918 


4895 


8.64 


2.01 


1.00 


1.06 


539 


247 


3 


1087 


173 





(1) Possible K secondary star 



Sp 



(R Q ) R CT (R ) R flt (R ) M fit (M ) T off 



MO.O 


0.543 


0.066 


0.490 


0.472 


M0.5 


0.528 


0.083 


0.488 


0.471 


M1.0 


0.429 


0.094 


0.480 


0.464 


Ml. 5 


0.443 


0.115 


0.465 


0.450 


M2.0 


0.468 


0.106 


0.445 


0.431 


M2.5 


0.422 


0.013 


0.420 


0.407 


M3.0 


0.415 


0.077 


0.391 


0.380 


M3.5 


0.361 


0.065 


0.359 


0.350 


M4.0 


0.342 


0.096 


0.326 


0.319 


M4.5 


0.265 


0.043 


0.292 


0.287 


M5.0 


0.261 


0.132 


0.258 


0.255 


M5.5 


0.193 


0.046 


0.226 


0.225 


M6.0 


0.228 


0.090 


0.195 


0.196 


M6.5 


0.120 


0.005 


0.168 


0.170 


M7.0 


0.178 


0.080 


0.145 


0.149 


M7.5 


0.118 


0.009 


0.126 


0.132 


M8.0 


0.137 


0.046 


0.114 


0.120 


M8.5 


0.110 


0.004 


0.109 


0.116 


M9.0 


0.108 


0.004 


0.112 


0.118 


M9.5 


0.111 


0.008 


0.124 


0.130 



above. In addition, it is difficult to assess the ag43 and metal- 

of the secondary stars in our WDMS sample. 

With the uncertainties on stellar parameters of sin- 

3843 gle M-dwarfs and the potential additional complications in 

3761 WDMS in mind, we decided to derive an "average" Sp — R 

3678 relation for M-dwarfs irrespective of their ages, metallicities, 

3596 anc j ac ti v ity levels. The primary purpose of this is to provide 

3514 distance estimates based on the flux scaling factors in Eq.lU 
3432 

^ a ^ so t° assess potential systematic peculiarities of the 

3267 seconc lary stars in the WDMS. 

3185 We ha ve compiled spectral typ e s and radii of field M- 
3103 dwarf s fr omlBerriman fc Rcid ( 198 fl).ICaillault fc Patterson! 



3090 (119901). ILeggett et al 



(2003) 
|200£ ) 



2938 Leto et al 
2856 
2773 
2691 
2609 
2527 
2444 



feOOCh. iLane etal 



(I 19991) 



Sceransan et al 

.Maceroni fc Montalbanl J2004t). ICreevev et al 
Pont et al l J200^. iRibasl d2006h. iBerger et al 



Table 3. Empirical Sp — R, Sp — M and Sp — T c ff relations (Rflt, 
Mfl t , T e fl ) found in this work. R mC an and R CT represent the mean 
radii and their standard deviation obtained from the sample of 
M-dwarfs described in Sect. 13.41 



2009 ). iBavless fc Orosj (|2006l ) and iBeattv et all (|2007h 
These data were separated into two groups, namely stars 
with directly measured radii (in eclipsing binaries or via 
2362 m terferometry) and stars with indirect radii determinations 
2281 ( e,e ' spectrophotometry) . We complemented this sample 
with spectral t ypes, masses, e f fectiv e temperatures and lu - 
minosit ie s from Dclfosse et al.l (119980 . ILeggett et al.l 1 2001). 
iBergerl (I2002TI. iGolimowski et all \2004 ). ICushing etaD 
(|2005h and iMontagnier et al l i2006h. calculating r adi i from 



47rflVr.ff 4 and/or ICaillault fc Pattersonl 's (ll99Cl ) 



mass-luminosity and mass-radius relations. 

Figure[7] shows our compilation of indirectly determined 



Besides the lack of extensive observational work on the 
Sp — R relation of single M-dwarfs, our need for an M-dwarf 
Sp — R relation in the context of WDMS faces a number of 
additional problems. A fraction of the WDMS in our sample 
have undergone a CE phase, and are now short-period bina- 
ries, in which the secondary star is tidally locked and hence 
rapidly rotating. This rapid rotation will enhance the stellar 
activity in a similar fashion to the short-period eclipsing M- 
dwarf binaries used in the M — R relation work mentioned 



In principle, an age estimate can be derived by adding the 
white dwarf cooling age to the main sequence life time of the 
white dwarf progenitor. This involves the use of an initia l mass - 
final mass relation for the white dwarf, e.g. iDobbie et al. 
which will not be strictly valid for those WDMS that underwent 
a CE evolution. Broadly judging from the distribution of white 
dwarf temperatures and masses in Fig. 1101 most WDMS in our 
sample should be older than 1 Gyr, but the data at hand does 
not warrant a more detailed analysis. 
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Figure 7. Top panel: indirectly measured radii of M-dwarfs vs 
spectral type. Our empirical Sp — R relation is given by a third 
order polynomial fit (solid line) to these data. Middle panel: 
mean radii and corresponding standard deviations obtained by 
averaging the radii in the top panel for each spectral type. Our 
Sp — R relation is again superimposed. Bottom panel: directly 
measured radii of M-dwarfs, again our empirical Sp — R rela- 
tion, the dashed li ne is the theoretical Sp — R relation from 
iBaraffe et all jl99Sh . M-dwarf radii from the eclipsing WDMS 
RRCae, NNSer, DECVn, RX J2130. 6+4710 and EC 13471-1258 
are shown as solid dots. 



radii as a function of spectral type (top panel) as well as 
those from direct measurements (bottom panel). A large 
scatter in radii is observed at all spectral types except for 
the very late M-dwarfs, where only few measurements are 
available. It is interesting that the amount of scatter is com- 
parable for both groups of M-dwarfs, those with directly 
measured radii and those with indirectly determined radii. 
This underlines that systematic effects intrinsic to the stars 
cause a large spread in the Sp — R relation even for the ob- 
jects with accurate measurements. In what follows, we use 
the indirectly measured radii as our primary sample, as it 
contains a larger number of stars and extends to later spec- 
tral types. The set of directly measured radii are used as 
a comparison to illustrate the Sp — R distribution of stars 
where the systematic errors in the determination of their 
radii is thought to be small. We determine an Sp — R rela- 
tion from fitting the indirectly determined radius data with 
third order polynomial, 

R = 0.48926 - 0.00683 Sp - 0.01709 Sp 2 + 0.00130 Sp 3 (2) 

The spectral type is not a physical quantity, and strictly 
speaking, this relation is only defined on the existing spectral 
classes. This fit agrees well with the average of the radii in 
each spectral class (Fig. [7] middle panel, where the errors 
are the standard deviation from the mean value). The radii 
from the polynomial fit are reported in Table|3j along with 
the average radii per spectral class. Both the radii from the 



polynomial fit and the average radii show a marginal upturn 
at the very latest spectral types, which should not be taken 
too seriously given the small number of data involved. 

We compare in Fig. [7] (bottom panel) the directly mea- 
sured radii with our Sp — R relation. It is apparent that also 
stars with well-determined radii show a substantial amount 
of scatter, and are broadly consistent with the empirical 
Sp — R relation determined from the indirectly measured 
radii. As a test, we included the directly measured radii 
in the fit described above, and did not find any significant 
change compared to the indirectly measured radii alone. 

For a final assessment on our empirical Sp — R 
relation, specifically in the context of WDMS, we have 
compiled from the literature the radii of M-dwa rfs in the 
eclipsing WDMS RR Cae dMaxted et all I2007T). NNSer 
dHaefner et all l2004h . DEC Vn (Ivan den Besselaar et ail 
l2007fi . RXJ21 30.6+4710 dMaxted et all |2004 ). and 
EC 13471-1258 dO'Donoghue et all l2003i r7Fig. bottom 
panel). Just as the accurate radii determined from inter- 
ferometric observations of M-dwarfs or from light curve 
analyses of eclipsing M-dwarf binaries, the radii of the 
secondary stars in WDMS display a substantial amount of 
scatter. 



3.4-1 Comparison with th e theoretical Sp-R relation from 
\Baraffe et all \l99&) 

We compare in the bottom panel of Fig. [7] our empirical 
Sp — R relation with the theoretical p r edicti on from the 
evolutionary sequences of IBaraffe et al.l (|l99ST ) , where the 
spectral type is based on the I — K colour of the PHOENIX 
stellar atmosphere models coupled to the stellar structure 
calculations. The theoretical Sp — R relation displays sub- 
stantially more curvature than our empirical relation, pre- 
dicting larger radii for spectral types <M2, and significantly 
smaller radii in the range M3-M6. The two relations con- 
verge at late spectral types (again, the upturn in the em- 
pirical relation for >M8.5 should be ignored as an artifact 
from our polynomial fit). The "kink" in the theoretical re- 
lation seen around M2 is thought to be a consequ ence of H2 
molecular dissociation dBaraffe fc Chabrier| [T996). The large 
scatter of the directly determined radii of field M-dwarfs as 
well of M-dwarfs in eclipsing WDMS could be related to 
two types of problem, that may have a common underlying 
cause. (1) In eclipsing binaries, the stars are forced to ex- 
tremely rapid rotation, which is thought to increase stellar 
activity that is likely to affect the stellar s tructure, gener- 
ally thought to lead to an increase in radius dSpruit fc Weisa 
Il98rj : iMullan fc MacDonaldll200ll ; IChabrier et alj2007l ). and 
(2) the spectral types in our compilations of radii are deter- 
mined from optical spectroscopy, and may differ to some 
extent from the spec tral type definition based on / — K 
colours as used in the IBaraffe et all dl99§t) models. Further- 
more, stellar activity is thought to affect not only the radii 
of the stars, but also their luminosity, surface temperatures, 
and hence spectral types. The effect of stellar activity is 
discussed in more detail in Sect . 14.71 

3.4-2 Sp — T c ff and Sp — M relations 

For completeness, we fitted the spectral type-mass data 
and the spectral type-effective temperature data compiled 
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from the literature listed above, and fitted the Sp — M and 
Sp — T c ft relations with a third-order polynomial and a first- 
order polynomial, respectively. The results from the fits are 
reported in Table|3] and will be used in this paper only for 
estimating upper limits to the orbital periods of our PCEBs 
(Sect. an d when discussing the possibility of stellar ac- 
tivity on the WDMS secondary stars in Sect . 14.71 



3.5 Distances 

The distances to the WDMS can be estimated from the best- 
fit flux scaling factors of the two spectral components. For 
the white dwarf, 

/wd / Rwd \ / \ 
~F~ = * \ ~r\ J (3) 

where / w d is the observed flux of the white dwarf, _F w d the 
astrophysical flux at the stellar surface as given by the model 
spectra, i? w a is the white dwarf radius and d w d is the dis- 
tance to the WD. For the secondary star, 



/sec 

p 

-t see 



/ R SC c \ " 
V dsec J 



(4) 



where / scc is the observed M-dwarf flux, F scc the flux 
at the stellar surface, and _R SCC and d ae c are the radius and 
the distance to the secondary respectively. 

The white dwarf radii are calculated from the best-fit 
T c ft and log<? as detailed in Sect. 13.31 The secondary star 
radii are taken from Table[3] for the best-fit spectral type. 
The uncertainties of the distances are based on the errors in 
R„d, which depend primarily on the error in logy, and in 
Rscc , where we assumed the standard deviation from Table[3] 
for the given spectral type. Table[S] lists the values d wd and 
d scc obtained for our PCEBs. The remaining 112 WDMS's 
distances can be found in the electronic edition (Tabled . 



4 DISCUSSION 

4.1 Ha vs Nai radial velocities 

As mentioned in Sect. [5] a few systems in Table[T]show con- 
siderable differences between their Ha and Nai radial veloc- 
ities. More specifically, while both lines clearly identify these 
systems as being radial velocity variable, and hence PCEBs 
or strong PCEB candidates, the actual radial velocities of 
Ha and Nai differ for a given SDSS spectrum by more than 
their errors. 

In close PCEBs with short orbital periods the Ha emis- 
sion is typically observed to arise from the hemisphere of 
the companion star facing the white dwarf. Irradiation from 
a hot white dwarf is the most plausible mechanism to ex- 
plain the anisotropic Ha emission, though also a number 
of PCEBs containing rather cool white dwarfs are known 
to exhibit concentrated Ha emissio n on the inner hemi- 
sphere of the com panion stars (e.g. Marsh fc Duckl 1 19961 ; 
iMaxted et ail feooei The anisotropy of the Ha emission re- 
sults in its radial velocity differing from other photospheric 
features that are (more) isotropically distributed over the 
companion stars, such as the Nai absorption. In general, 
the Ha emission line radial velocity curve will then have a 
lower amplitude than that of the Nai absorption lines, as 



Ha originates closer to the centre of mass of the binary sys- 
tem. In addition, the strength of Ha can vary greatly due to 
different geometric projections in high inclination systems. 
More complications are added in the context of SDSS spec- 
troscopy, where the individual spectra have typical exposure 
times of 45-60min, which will result in the smearing of the 
spectral features in the short-period PCEBs due to the sam- 
pling of different orbital phases. This problem is exacerbated 
in the case that the SDSS spectrum is combined from ex- 
posures taken on different nights (see Sect .[2)1 . Finally, the 
Ha emission from the companion may substantially increase 
during a flare, which will further enhance the anisotropic na- 
ture of the emission. 

Systems in which the Ha and Nai ra- 
dial velocities differ by more than 2 a are: 
SDSS J005245.11-005337.2, SDSS J024642.55+004137.2, 
SDSSJ030904.82-010100.8, SDSS J031404. 98-011136. 6, and 
SDSS J172406. 14+562003.0. Of these, SDSS J0246+0041, 
SDSS J0314-0111, and SDSS J1724+5620 show large- 
amplitude radial velocity variations and substantial changes 
in the equivalent width of the Ha emission line, suggesting 
that they are rather short orbital period PCEBs with 
moderately high inclinations, which most likely explains 
the observed differences between the observed Ha and Nai 
radial velocities. Irradiation is also certainly important 
in SDSS J1724+5620 which contains a hot (~ 36000 K) 
white dwarf. SDSS J0052-0053 displays only a moderate 
radial velocity amplitude, and while the Ha and Nai radial 
velocities display a homogeneous pattern of variation (Fig. [3] 
and 2| , Ha appears to have a larger amplitude which is 
not readily explained. Similar discrepancies have been 
observed e.g. in the close magnetic WDMS binary WXLMi, 
and were thought to be related to a ti me- variable chang e 
in the location of the Ha emission (|Vogel et al.l I2007T ). 
Finally, SDSS J0309-0101 is rather faint (g = 20.4), but has 
a strong Ha emission that allows reliable radial velocity 
measurements that identify the system as a PCEB. The 
radial velocities from the Nai doublet are more affected by 
noise, which probably explains the observed radial velocity 
discrepancy in one out of its seven SDSS spectra. 

4.2 Upper limits to the orbital periods 

The radial velocities of the secondary stars follow from Ke- 
pler's 3rd law and depend on the stellar masses, the orbital 
period, and are subject to geometric foreshortening by a fac- 
tor sini, with i the binary inclination with regards to the 
line of sight: 



(M wd sinti) 3 
(M wd + M soc ) 2 



P uK^ 

1 or b ± v st 

2ttG 



(5) 



with .Ksec the radial velocity amplitude of the secondary 
star, and G the gravitational constant. This can be rear- 
ranged to solve for the orbital period, 

27rG(M wd sini) 3 



(M wd + M sec ) 2 KL 



(6) 



From this equation, it is clear that assuming i — 90° gives 
an upper limit to the orbital period. 

The radial velocity measurements of our PCEBs and 
PCEB candidates (Table[TJ sample the motion of their com- 
panion stars at random orbital phases. However, if we as- 
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Figure 8. Comparison of the white dwarf effective temperatures, distances base d on the white dwarf fi t, and the spectral types of the 
secondary stars determined from our fits (Sect, l3~2l I3,3l and Tabled, and those of lRavmond et al I J2003h . Top panels, from left to right: 
the ratio in T ef j , the ratio in d, and the difference in the secondary's spectral types from the two studies as a function of the white dwarf 
temperature. 



sume that the maximum and minimum values of the ob- 
served radial velocities sample the quadrature phases, e.g. 
the instants of maximum radial velocity, we obtain lower 
limits to the true radial velocity amplitudes of the compan- 
ion stars in our systems. From Eq.HJ] a lower limit to -Ksec 
turns into an upper limit to P or b- 

Hence, combining the radial velocity information from 
Table[T]with the stellar parameters from TableJSJ we deter- 
mined upper limits to the orbital periods of all PCEBs and 
PCEB candidates, which range between 0.46-7880 d. The 
actual periods are likely to be substantially shorter, espe- 
cially for those systems where only two SDSS spectra are 
available and the phase sampling is correspondingly poor. 
More stringent constraints could be obtained from a more 
complex exercise where the mid-exposure times are taken 
into account - however, given the fact that many of the 
SDSS spectra are combined from data taken on different 
nights, we refrained from this approach. 

4.3 The fraction of PCEB among the SDSS 
WDMS binaries 

We have measured the radial velocities of 101 WDMS which 
have multiple SDSS spectra, and find that 15 of them clearly 
show radial velocity variations, three additional WDMS are 
good candidates for radial velocity variations (see Tabled . 
Taking the upper limits to the orbital periods at face value, 
and assuming that systems with a pe riod < 300 d have un- 
dergone a CE (|Willems fc Kolbll2004 see also Sect.© 17 of 
the systems in TableQ] qualify as PCEBs, implying a PCEB 
fraction of ~15% in our WDMS sample, which is in rough 
agreement with the predictions by the population model 



of IWillems fc Kolbl l|2004h . However, our value is likely to 
be a lower limit on the true fraction of PCEBs among the 
SDSS WDMS binares for the following reasons. (1) In most 
cases only two spectra are available, with a non-negligible 
chance of sampling similar orbital phases in both observa- 
tions. (2) The relatively low spectral resolution of the SDSS 
spectroscopy (A/AA ~ 1800) plus the uncertainty in the flux 
calibration limit the detection of significant radial velocity 
changes to ~ 15kms _1 for the best spectra. (3) In binaries 
with extremely short orbital periods the long exposures will 
smear the Nai doublet beyond recognition. (4) A substan- 
tial number of the SDSS spectra are combined, averaging 
different orbital phases and reducing the sensitivity to radial 
velocity changes. Follow-up observations of a representative 
sample of SDSS WDMS with higher spectral resolution and 
a better defined cadence will be necessary for an accurate 
determination of the fraction of PCEBs. 

4.4 Comparison with iRavmond et al.l {2003) 

In a previous study, IRavmond et al.l l|2003h determined 
white dwarf temperatures, distance estimates based on the 
white dwarf fits, and spectral types of the companion star 
for 109 SDSS WDMS. They restricted their white dwarf fits 
to a single gravity, log g = 8.0, and a white dwarf radius of 
8 x 10 8 cm (corresponding to M w d = 0.6 M©), which is a fair 
match for the majority of systems (see Sect. l4.6l below'). Our 
sample of WDMS with two or more SDSS spectra has 28 ob- 
jects in common with Raymond's list, sufficient to allow for 
a quantitative comparison between the two different meth- 
ods used to fit the data. As we fitted two or more spectra for 
each WDMS, we averaged for this purpose the parameters 
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Figure 9. Comparison of the white dwarf effective temperature s and surface graviti es and the spectral types of the secondary stars 
determined from our fits (Sect.[32]EI]and Table|5j, and those of lSilvestri et al.|j2006h . Top panels, from left to right: the WD effective 
temperature and surface gravity ratios, and the difference in the secondary's spectral types from the two studies as a function of the 
white dwarf temperature. 



obtained from the fits to individual spectra of a given ob- 
ject, and propagated their errors according ly. We find that 
~ 2/3 of the temperatures determined by iRavmond et al.l 
(2003) agree with ours at the ~ 20 per cent level, with 
the remaining being different by up to a factor two (Fig.[S] 
left panels). This fairly large disagreement is most likely 
caused by the simplified fitting Raymond et al. adopted, 
i.e. fitting the white dwarf models in the wavelength range 
3800-5000 A, neglecting the contribution of the compan- 
ion star. The spectral types of the companion stars from 
our work and lRavmond et alj l|2003l ) agree mostly to within 
±1.5 spectral classes, which is satisfying given the compos- 
ite nature of the WDMS spectra and the problems associ- 
ated with their spectral decomposition (Fig.|H] right panels). 
The biggest discrepancy shows up in the distances, with the 
Raymond et al. distances being systematically lower than 
ours (Fig. [5] middle panels). The average of the factor by 
which Raymond et al. underpredict the distances is 6.5, 
which is close to 2-7T, suggesting that the authors may have 
misinterpreted the flux definition of the m odel atmosphere 
code they used (TLUSTY/SYNSPEC from lHubenv fc Land 
1 19951 . which outputs Eddington fluxes) , and hence may have 
used a wrong constant in the flux normalisation (Eq.[3J). 

4.5 Comparison with ISilvestri et all (| 20061 ) 

Having developed an independent method of determining 
the stellar parameters for WDMS fro m their SPSS s pectra , 
we compared our results to those of ISilvestri et al] |2006). 
As in Sect. 14.41 above, we average the parameters obtained 
from the fits to the individual SDSS spectra of a given ob- 
ject. Figure|9]shows the comparison between the white dwarf 



effective temperatures, surface gravities, and spectral types 
of the secondary stars from the two studies. Both studies 
agree in broad terms for all three fit parameters (Fig. [9] 
bottom panels). Inspecting the discrepancies between the 
two independent sets of stellar parameters, it became evi- 
dent that relatively large disagreements are most noticeably 
found for T g < 20 000 K, with differences in T e g of up to 
a factor two, an order of magnitude in surface gravity, and 
a typical difference in spectral type of the secondary of ±2 
spectral classes. For higher temperatures the differences be- 
come small, with nearly identical values for T c g, logg agree- 
ing within ±0.2 magnitude, and spectral types differing by 
±1 spectral classes at most (Fig. [9] top panels). We inter- 
pret this strong disagreement at low to intermediate white 
dwarf temperatures to the ambiguity between hot and cold 
solutions described in Sect. 13.31 

A quantitative judgement of the fits in ISilvestri et alj 

(2006) is difficult, as the authors do not provide much de- 
tail on the method used to decompose the WDMS spectra, 
except for a single example in their Fig. 1. It is worth noting 
that the M dwarf component in that figure displays constant 
flux at A < 6000 A, which seems rather unrealistic for the 
claime d spectral type of M5. Unfortunately, ISilvestri et al.l 
(2006) do not list distances implied by their fits to the white 
dwarf and main sequence components in their WDMS sam- 
ple, which would provide a test of internal consistency (see 

Sect.gTJ). 

We also investigated the svstems lSilvestri et alj 's (2006) 
method failed to fit, and found that we were able to deter- 
mine reasonable parameters for the majority of them. It ap- 
pears that our method is more robust in cases of low signal- 
to-noise ratio, and in cases where one of the stellar compo- 
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Figure 10. White dwarf mass, Sp types of the secondaries, effec- 
tive temperature and logg histograms obtained from the SDSS 
WDSS sample. Excluded are those systems with individual WD 
masses, T e g, and log g associated to relative errors larger than 25 
per cent. 



nents contributes relatively little to the total flux. Examples 
of the latter are SDSS J204431. 45-061440. 2, where an MO 
secondary star dominates the SDSS spectrum at A > 4600 A, 
or SDSS J172406. 14+562003.1, which is a close PCEB con- 
taining a hot white dwarf and a low-mass companion. An in- 
dependent analysis of the entire WDMS sample from SDSS 
appears therefore a worthwhile exercise, which we will pur- 
sue elsewhere. 



4.6 Distribution of the stellar parameters 

Having determined stellar parameters for each individual 
system in Sect.|3j we are looking here at their global distri- 
bution within our sample of WDMS. Figure. [101 shows his- 
tograms of the white dwarf effective temperatures, masses, 
logg, and the spectral types of the main-sequence compan- 
ions. 

As in Sect. 14.41 and 14.51 above, we use here the aver- 
age of the fit parameters obtained from the different SDSS 
spectra of each object. Furthermore, we exclude all sys- 
tems with relative errors in their white dwarf parameters 
{T w d, log g, M w d) exceeding 25 per cent to prevent smear- 
ing of the histograms due to poor quality data and/or fits, 
which results in 95, 81, 94, and 38 WDMS in the his- 
tograms for the companion spectral type, logg, T wd , and 
M w( j, respectivel y. In broad terms, our resu lts are consis- 
tent w ith those of lRavmond et all (|2003t ) and lSilvestri et all 
(2006): the most frequent white dwarf temperatures are be- 
tween 10 000-20 000 K, white dwarf masses cluster around 
M w( j — O.6M0, and the companion stars have most typi- 
cally a spectral type M3-4, with spectral types later than 
M7 or earlier than Ml being very rare. 



At closer inspection, the distribution of white dwarf 
masses in our sample has a more pronounce d tail towards 
lower masses compared to the distribution in lSilvestri et al.l 
(2006). A tail of lower- mass white dwarfs, peaking around 
0.4 Mq is observe d also in well-studie d samples of single 
white dwarfs (e.g. iLiebert et al.l 12005). and is interpreted 
as He-core whi te dwarfs descendi ng from evolution in a bi- 
nary star (e.g. iMarsh et all 1 19951 ). In a sample of WDMS, 
a significant fraction of systems will have undergone a CE 
phase, and hence the fraction of He-core white dwarfs among 
WDMS is expected to be larger than in a sample of single 
white dwarfs. 

Also worth noting is that our distribution of com- 
panion star spectral types is relatively flat between M2- 
M4, more similar to the d istribution of single M-dwarfs 
in SDSS dWest et alj 120041 ') than the companion stars in 
ISilvestri et al. 1 20061 ). More generally speaking, the cut-off 
at early spectral types is due to the fact that WDMS with 
K-type companions can only be identified from their spec- 
tra/colours if the white dwarf is very hot - and hence, very 
young, and correspondingly only few of such systems are in 
the total SDSS WDMS sample. The cut-off seen for low- 
mass companions is not so trivial to interpret. Obviously, 
very late-type stars are dim and will be harder to be de- 
tected against a moderately hot white dwar f, such a bias 
was discussed bv lSchreiber fc Gansickel (|2003l ) for a sample 
of ~ 30 well-studied WDMS which predominantly originated 
from blue-colour ( = hot white dwarf) surveys. However, old 
WD MS with cool white dwarfs sh ould be much more com- 
mon ( Sch reiber fc Gansickel [2003). and SDSS, sampling a 
much broader colour space than previous surveys, should be 
able to identify WDMS containing cool white dwarfs plus 
very late type companions. The relatively low frequency of 
such systems in the SDSS spectroscopic data base suggests 
that either SDSS is not efficiently targeting those systems 
for spectroscopic follow-up, or that they are rare in the first 
place, or a combination of both. A detailed disc ussion is be- 
yond the scope of this paper, but we note that iFarihi et alj 
(|2005h have constructed the relative distribution of spec- 
tral types in the local M/L dwarf distribution, which peaks 
around M3-4, and steeply declines towards later spectral 
types, suggesting that late- type companions to white dwarfs 
are intrinsically rare. This i s supported independently by 
iGrether fc Lineweaverl (|2006), who analysed the mass func- 
tion of companions to solar-like stars, and found that it 
steeply decreases towards the late end of the main sequence 
(but rises again for planet-mass companions, resulting in the 
term "brown dwarf desert"). 

An assessment of the stellar parameters of all WDMS 
in SDSS DR5 using our spectral decomposition and white 
dwarf fitting method will improve the statistics of the dis- 
tributions presented here, and will be presented in a future 
paper. 



4.7 Stellar activity on the secondary stars? 

As outlined in Sect . 13.51 the scaling factors used in the mod- 
elling of the two spectral components of each WDMS provide 
two independent estimates of the distance to the system. In 
principle, both estimates should agree within their errors. 
Figure[TT] compares the white dwarf and secondary star dis- 
tance estimates obtained in Sect. 13.51 where the distances 
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Figure 11. Comparison of d Bec and d w< j obtained from our spectral decomposition and white dwarf fits to the SDSS spectra. Approx- 
imately a third of the systems have d Bec ^ rf w d- The left panel splits the sample according to the mass of the white dwarfs, while the 
right panel divides the sample according to the spectral types of the secondaries. In both panels systems that we identify as PCEBs from 
radial velocity variations in their SDSS spectra are shown in red. 




Figure 12. Left panel: the di s tance s implied by the spectral decomposition were calculated by using the Sp — R relation predicted 
by the models of Baraffe ct al. (1998:), instead of our empirical Sp — R relation. Right panel: the spectral types of the secondary stars 
were adjusted by 1-2 spectral classes to achieve d wd = d scc . Only three systems can not be reconciled in this way, and are discussed 
individually in the text. We suggest that stellar activity in some WDMS may change the spectral type of their secondary stars, being 
equivalent to a change in surface temperature by a few 100 K. 



obtained from the individual SDSS spectra of a given ob- 
ject were averaged, and the errors accordingly propagated. 
In this plot, we exclude systems with relative errors in ci w d 
larger than 25 per cent to avoid cluttering by poor S /N data. 
The relative error in d sec is dominated by the scatter in the 
Sp — R relation, which represents an intrinsic uncertainty 
rather than a statistical error in the fit, and we therefore 
did not apply any cut in d scc . Taking the distribution of 
distances at face value, it appears that about 2/3 of the 
systems have d scc ~ d w( j within their la errors, as expected 
from purely statistical errors. However, there is a clear trend 
for outliers where d scc > d w d- We will discuss the possible 
causes and implications in the following sections. 



4-7.1 Possible causes for d sec 7^ d w d 

We identify a number of possible causes for the discrepancy 
between the two independent distance estimates observed in 
~ 1/3 of the WDMS analysed here. 

(1) A tendency for systematic problems in the white 
dwarf fits? d scc > d w d could be a result of too small white 
dwarf radii for a number of systems, i.e. too high white dwarf 
masses. We therefore identify in the left panel of Fig. 1111 
those systems with massive (> 0.75 M©) white dwarfs. It is 
apparent that the outliers from the d S ec = d w d relation do 
not contain a large number of very massive white dwarfs. 

(2) Problems in determining the correct spectral type 
of the secondary? If the error on the spectral type of the 
companion star determined from the spectral decomposition 
is larger than ±0.5, as assumed in Sect. 13.21 a substantial 
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deviation from d scc = d wd would result. However, as long 
as this error is symmetric around the true spectral type, it 
would cause scatter on both sides of the d sec ~ dwd relation. 
Only if the determined spectral types were consistently too 
early for ~ 1/3 of the systems, the observed preference for 
outliers at d BCC > d wd could be explained (see Sect. 14.7.21 
below for a hypothetical systematic reason for spectral types 
that are consistently too early). 

(3) Problems in the spectral type-radius relation? As dis- 
cussed in Sect. 13.41 the spectral type-radius relation of late 
type stars is not particularly well defined. The large scat- 
ter of observed radii at a given spectral type is taken into 
account in the errors in d sec . If those errors were underes- 
timated, they should cause an approximatively symmetric 
scatter of systems around d soc = d w d, which is not observed 
(the Sp — R relation being non-linear lead to asymmetric 
error bars in the radius for a given symmetric error in the 
spectral type, however, over a reasonably small range in the 
spectral type this effect is negligible). A systematic problem 
over a small range of spectral types would result in a con- 
centration of the affected spectral types among the outliers. 
For this purpose, we divide our sample into three groups of 
secondary star spectral classes, Sp > 5, 3 < Sp < 5, and 
Sp < 3 (Fig. 1111 right panel). The outliers show a slight 
concentration towards early types (Sp < 3) compared to 
the distribution of secondary star spectral types in the total 
sample (Fig. [TTJJ). 

To explore the idea that our empirical Sp — R relation 
is simply inadequate, we calculated a new set of secondary 
star distances , usin g the theoretical Sp — R relation from 
iBaraffe et al.l |l998;) (see Fig.0 bottom panel), which are 
shown in the left panel of Fig. 1121 The theoretical Sp — R 
relation implies smaller radii in the range M3-M6, but the 
difference with our empirical relation is not sufficient enough 
to shift the outlying WDMS onto the d scc — d w d relation. 
For spectral types earlier than M2.5, our empirical Sp — 
R relation act ually gives smaller radii than the theoretical 
Baraffe et al. Il998l relation, so that using the theoretical 
Sp — R actually exacerbates the d scc > d w( j problem. 

(4) A relationship with close binarity? The fraction of 
PCEBs among the outliers is similar to the fraction among 
the total sample of WDMS (Fig. llljl . hence it does not ap- 
pear that close binarity is a decisive issue. 

(5) An age effect? Late type stars take a long time to 
contract to their zero age main sequence (ZAMS) radii, and 
if some of the WDMS in our sample were relatively young 
objects, their M-dwarfs would tend to have larger radii than 
ZAMS radii. As briefly discussed in Footnote^ the major- 
ity of the WDMS in our sample are likely to be older than 
~ 1 Gyr, and the outliers in Fig. lllll2l do not show any pref- 
erence for hot or massive white dwarfs, which would imply 
short cooling ages and main sequence life times. 

4-7.2 Could stellar activity affect Sp scc ? 

None of the points discussed in the previous section conclu- 
sively explains the preference for outliers having d scc > d w( j. 
If we assume that the problem rests in the determined prop- 
erties of the secondary star, rather than those of the white 
dwarf, the immediate implication of d SC c > d wd is that the 
assumed radii of the secondary stars are too large. As men- 
tioned above and shown in Fig. 1121 this statement does not 



strongly depend on which Sp — R relation we use to deter- 
mine the radii, eit her ou r empirical relation or the theoret- 
ical Baraffe et al. (I1998T ) relation. Rather than blaming the 
radii, we explore here whether the secondary star spectral 
types determined from our decomposition of the SDSS spec- 
tra might be consistently too early in the outlying systems. 
If this was the case, we would pick a radius from our Sp — R 
relation that is larger than the true radius of the secondary 
star, resulting in too large a distance. In other words, the 
question is: is there a mechanism that could cause the spec- 
tral type of an M-star, as derived from low-resolution optical 
spectroscopy, to appear too early? 

The reaction of stars to stellar activity on their sur- 
face, also referred to as spottedness is a complex phe- 
nomeno n that is not fully understood. Theoretical stud- 



Spruit fc Weiss! 1 19861 ; iMullan fc MacDonaldl l200ll ; 
IChabrier et all 2007) agree broadly on the following points: 
(1) the effect of stellar activity is relatively weak at the 
low-mass end of the main sequence (M < 0.3 M0), where 
stars are conve ntionally thought to become fully convec- 
tiye (though, see lMullan fc MacDonaldll200ll ; IChabrier et all 
2007 for discussions on how magnetic fields may change that 
mass boundary), (2) stellar activity will result in an increase 
in radius, and (3) the effective temperature of an active star 
is lower than that of an unspotted star. 

Here, we briefly discuss the possible effects of stellar ac- 
tivity on the spectral type of a star. For this purpose, it is 
important not to confuse the effective temperature, which is 
purely a definition coupled to the luminosity and the stellar 
radius via L = 4irR 2 aT^ g (and hence is a global property of 
the star), and the local temperature of a given part of the 
stellar surface, which will vary from spotted areas to inter- 
spot areas. In an unspotted star effective and local temper- 
ature are the same, and both colour and spectral type are 
well-defined. As a simple example to illustrate the difference 
between effective temperature and colour in an active star, 
we assume that a large fraction of the star is covered by 
zero-temperature, i.e. black spots, and that the inter-spot 
temperatu re is the same a s that of the unspotted star. As 
shown by Cha brier et al.1 (|2007l ). assuming constant lumi- 
nosity requires the radius of the star to increase, and the 
effective temperature to drop. Thus, while intuition would 
suggest that a lower effective temperature would result in a 
redder colour, this ficticious star has exactly the same colour 
and spectral type as its unspotted equivalent - as the black 
spots contribute no flux at all, and the inter-spot regions 
with the same spectral shape as the unspotted star. 

Obviously, the situation in a real star will be more com- 
plicated, as the spots will not be black, but have a finite tem- 
perature, and the star will hence have a complicated temper- 
ature distribution over its surface. Thus, the spectral energy 
distribution of such a spotted star will be the superposition 
of contributions of different temperatures, weighted by their 
respective covering fraction of the stellar surface. Strictly 
speaking, such a star has no longer a well-defined spec- 
tral type or colour, as these proper ties will depend on the 
wavelength range that is observed. ISpruit fc Weissl (Il986l ) 
assessed the effect of long-term spottedness on the temper- 
ature distribution on active stars, and found that for stars 
with masses in the range 0.3 — 0.6 Mq the long-term effect of 
spots is to increase the temperature of the inter-spot regions 
by ~ 100 — 200 K (compared to the effective temperature of 
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the equivalent unspotted star), wheras the inter-spot tem- 
perature of spotted l ower-mass stars remains unchanged. 
ISpruit fc Weissl (| 19861 } also estimated the effects of stellar 
activity on the colours of stars, but given their use of sim- 
ple blackbody spectra, these estimates are of limited value. 
As a general tendency, the hotter (unspotted) parts of the 
star will predominantly contribute in the blue end of the the 
spectral energy distribution, the cooler (spotted) ones in its 
red end. As we determine the spectral types of the secondary 
stars in the SDSS WDM S from optic al ( = blue) spectra, and 
taking the results of lSpruit fc Weissl at face value, it appears 
hence possible that they are too early compared to unspot- 
ted stars of the same mass. A full theoretical treatment of 
this problem would involve calculating the detailed surface 
structure of active stars as well as appropriate spectral mod- 
els for each surface element in order to compute the spatially 
integreated spectrum as it would be observed. This is clearly 
a challenging task. 

Given that theoretical models on the effect of stellar 
activity have not yet converged, and are far from making 
detailed predictions on the spectroscopic appearence of ac- 
tive stars, we pursue here an empirical approach. We as- 
sume that the discrepancy d scc > d w( j results from picking 
a spectral type too early, i.e. we assume that the secondary 
star appears hotter in the optical spectrum that it should 
for its given mass. Then, we check by how much we have 
to adjust the spectral type (and the corresponding radius) 
to achieve d scc = d w d within the errors. We find that the 
majority of systems need a change of 1-2 spectral classes, 
which corresponds to changes in the effective temperature 
of a few hundred degrees only, in line with the calculations 
of ISpruit fc Weissl (|l986l ). Bearing in mind that what we 
see in the optical is the surface temperature, and not the 
effective temperature, compari ng this to the surface tem- 
perature changes calculated bv lSpruit fc We iss (1986), and 
taking into account that we ignored in this simple approach 
the change in radius caused by a large spottedness, it ap- 
pears plausible that the large deviations from d scc — d wt j 
may be related to stellar activity on the secondary stars. 

There are three WDMS where a change of more than 
two spectral classes would be necessary: SDSSJ032510.84- 
011114.1, SDSSJ093506.92+441107.0, and SDSSJ143947.62- 
010606.9. SDSSJ143947.62-010606.9 contains a very hot 
white dwarf, and the secondary star may be heated if 
this system is a PCEB. Its two SDSS spectra reveal no 
significant radial velocity variation, but as discussed in 
Sect. [5] the SDSS spectra can not exclude a PCEB nature 
because of random phase sampling, low inclination and 
limited spectral resolution. SDSSJ032510.84-011114.1 and 
SDSSJ093506.92+441107.0 could be short-period PCEBs, 
as they both have poorly define Nai absorption doublets, 
possibly smeared by orbital motion over the SDSS expo- 
sure (see Sect. [2}. In a close binary, their moderate white 
dwarf temperatures would be sufficient to cause noticeable 
heating of the secondary star. We conclude that our study 
suggests some anomalies in the properties of ~ 1/3 of the M- 
dwarf companions within the WDMS sample analysed here. 
This is in line with previous detailed studies reporting the 
anomalous behaviour of the main seq uence compan i ons in 
P CEBs and cat a clysm ic variables, e.g. [o 'Brien et all (120011 ) 
or lNavlor et~ai1 (|2005h . 



4.8 Selection effects among the SDSS WDMS. 

Selection effects among the WDMS found by SDSS with re- 
spect to the spectral type of their main-sequence component 
can be deduced from the right panel of Fig. 1111 No binaries 
with econdary spectral types later than M5 are found at 
distances larger than ~ 500 pc. Because of their intrinsic 
faintness, such late- type secondary stars can only be seen 
against relatively cool white dwarfs, and hence the large 
absolute magnitude of such WDMS limits their detection 
within the SDSS magnitude limit to a relatively short dis- 
tance. Hot white dwarfs in SDSS can be detected to larger 
distances, and may have undetected late-type companions. 
There are also very few WDMS with secondary stars ear- 
lier than M3 within 500 pc. In those systems, the secondary 
star is so bright that it saturates the z, and possibly the i 
band, disqualifying the systems for spectroscopic follow up 
by SDSS. While these selection effects may seem disheart- 
ing at first, it will be possible to quantitatively correct them 
based on predicted colours of WDMS binaries and the in- 
formation available within the SDSS project regarding pho- 
tometric properties and spectroscopic selection algorithms. 



5 CONCLUSIONS 

We have identified 18 PCEBs and PCEB candidates among 
a sample of 101 WDMS for which repeat SDSS spectroscopic 
observations are available in DR5. From the SDSS spectra, 
we determine the spectral types of the main sequence com- 
panions, the effective temperatures, surface gravities, and 
masses of the white dwarfs, as well as distance estimates 
to the systems based both on the properties of the white 
dwarfs and of the main sequence stars. In about 1/3 of the 
WDMS studied here the SDSS spectra suggest that the sec- 
ondary stars have either radii that are substantially larger 
than those of single M-dwarfs, or spectral types that are too 
early for their masses. Follow-up observations of the PCEBs 
and PCEB candidates is encouraged in order to determine 
their orbital periods as well as more detailed system param- 
eters. Given the fact that we have analysed here only ~ 10 
per cent of the WDMS in DR5, it is clear SDSS holds the 
potential to dramatically improve our understanding of CE 
evolution. 
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Table 5: WD masses, effective temperatures, surface gravities, spectral 
types and distances of the SDSS PCEBs identified in Sect.O as deter- 
mined from spectral modelling. The stellar parameters for the remaining 
112 WDMS binaries can be found in the electronic edition of the paper. 
We quot e by and s and e tho se sys tems which have been studied previ- 
ously bv lSilvestri et all l|2006l ) and lEisenstein et al.l (|2006l ) . repectively. 
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106 


6 


639 


288 




224139.02+002710.9 


53261 


1901 


471 


12681 


495 


8.05 


0.15 


0.64 


0.09 


369 


35 


4 


381 


113 


e 




52201 


674 


625 


13745 


1644 


7.66 


0.36 


0.43 


0.19 


524 


108 


4 


378 


112 




233928.35-002040.0 


53357 


1903 


264 


15071 


1858 


8.69 


0.33 


1.04 


0.18 


416 


112 


4 


530 


157 


e 




52525 


682 


159 


12536 


2530 


7.92 


0.79 


0.56 


0.48 


655 


291 


4 


528 


157 




234534.49-001453.7 


52524 


683 


166 


19193 


1484 


7.79 


0.31 


0.51 


0.17 


713 


132 


4 


1058 


314 


s,e 



Continued on Next Page. 
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Table 5 - Continued 



SDSS J 


MJD 


plate 


fiber 


Teff(K) 


err log g err 


M(Mq) err 


dwd(pc) err Sp 


4oc(pc) 


err flag notes 


235020.76-002339.9 


53357 
51788 
52523 


1903 

386 
684 


103 
228 
226 


18974 


730 7.98 0.15 


0.61 0.09 


652 62 4 

5 
5 


1155 
504 
438 


343 

254 6 
22 



(1) T c ff less than 6000; (2) Noisy spectra; (3) Cold WD; (4) Diffuse background galaxy in the SDSS image; (5) Reflection 
effect; (6) Some blue excess, WD? 



